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Abstract 

Ultra-cold neutrons (UCN) are neutrons with energies so low they can be stored in material bot- 
tles and magnetic traps. They have been used to provide the currently most accurate experiments 
on the neutron life time and electric dipole moment. UCN can be produced in superfluid Helium 
at significantly higher densities than by other methods. The predominant production process is 
usually by one phonon emission which can only occur at a single incident neutron energy because 
of momentum and energy conservation. However UCN can also be produced by multiphonon pro- 
cesses. It is the purpose of this work to examine this multiphonon production of UCN. We look 
at several different incident neutron spectra, including cases where the multiphonon production is 
significant, and see how the relative importance of multiphonon production is influenced by the 
incident spectrum. 

PACS numbers: 78.70.Nx, 67.20+k, 61.12.Ex, 29.25.Dz 
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Introduction 

The production of Ultra-Cold Neutrons (UCN)0 in superfluid Helium and their subse- 
quent storage allowing the build up of densities much greater than would be achievable 
at thermal equilibrium with the Helium (superthermal UCN source)! is the foundation for 
several current and proposed experiments with UCNiffli. Although calculations of the pro- 
duction rate have been presented previouslyifl'0 the early calculations concentrated on the 
production due to the resonant (8.9A) phonons. Measurements! have shown a rather nar- 
row peak in the production cross section as expected. In the course of the first calculation! 
estimates had been made of the contribution due to higher energy incident neutrons down- 
scattering on the multiphonon excitation tail in 5* (g, u) , and these neutrons were found 
to contribute about the same amount of UCN production as the resonant neutrons for a 
reasonable assumed spectrum. However as the initial work was part of a proposal for an 
experiment it was felt that it would be more conservative not to include this in the initial 
estimates of UCN density i.e. to keep this in hand as a factor of safety for the predicted 
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intensity for the new source. However the fact that multiphonon processes would increase 
the UCN production rate was mentioned. 

In recent years the comparison between theory and experiment has improved signifi- 
cantly and there are several proposals to operate a superthermal Helium UCN source with a 
monochromatic beam of 8.9 A neutrons, so the question of the production rate due to multi- 
phonon processes has some practical importance. In addition, due to a vigorous experimental 
program!, much more detailed information concerning S (g, cu) has become available. For 
these reasons it has become necessary to re-examine the situation. W. Schott, in particular, 
has emphasized the importance of the UCN production due to higher energy neutrons. 

This work is concerned with calculations of the UCN production rate when a volume of 
superfluid He'' at temperatures of < 0.57^ is exposed to a neutron flux. The steady state 
UCN density in the source is given by 

PUCN = Pt (1) 

where P is the production rate, UCN/cm^ / sec and r is the UCN lifetime against all losses 
in the Helium filled storage chamber. Lifetimes > 500 seconds should be achievable. The 
storage time in superfluid Helium has been measured 

I. PRODUCTION OF UCN 
A. Introduction 

A neutron at rest can absorb energy hw and momentum hq with 

u = hq^/2m = (2) 

with {a = 4.14 mev/A^^) and contrarily a neutron with this energy and momentum can 
come to rest after transferring its energy and momentum to the superfluid. For single phonon 
interactions, which are usually dominant, the superfluid can only exchange quantities of 
energy and momentum that are related by the dispersion curve 

UJ = UJ (q) = cq (3) 

where the last relation is an approximation to simplify the discussion. Then neutrons can 
only come to rest by emission of a single phonon if they have the resonant energy E* given 
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by the intersection of equations (0) and (|^) 

uj (q) = cq = hq^ /2m 
q* = 2mc/h 

In addition to one phonon interactions there are multiple phonon interactions occurring at 
higher input energies that can also produce UCN. These are the main subject of this work. 

UCN production in Helium can be calculated as follows: The differential cross section for 
neutron scattering is given by the Fourier transform of the van Hove correlation function 
S{q,u;)E: 



^ = b''^S{q,u)dn 



and S{q,uj) has been measured in great detail!. 



(4) 



dVL = 27ismede = 27i 



qdq 
kik2 



Then 



Since the limits on q are: 

ki - k2 < q < ki + k2 
k2 = ku ki q ki {hku is the UCN momentum) 

we have dq = 2ku. Then 



(5) 



(6) 



where we assumed S{q,u) is constant over the narrow range dq. 



_ ^(fef fej) _ a (1.2 _ t.2\ ^ ai.2 



The UCN production rate is given by 



(7) 



P{Eu)dEu 



'-^Nj-^iE,-^E.)dE, 
dh duo 



dEu 



(8) 
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then 



dE 
d^{k 



—S{ki, uj 



dE 



—S{ki, UJ 



akj 
ak\ 



)dki 
)dki 



k^dku 



UCN/cmy sec (9) 



where E^, k^ are the critical UCN energy and wave vector of the walls of the storage chamber 



and 



d'I'(fci) 



dE dE 

Note S{ki,uj) is in {mev)~^ 



is the energy spectrum of the incoming flux. 



B. Single phonon production of UCN 



The one phonon production ratea'Bllj is found by evaluating the integral in equ.(|^) over 
the one phonon peak. (Note that k* = .7 A^^) . 

Pi„ = NA7rb^S*a(3^ UCN/cmV sec (10) 

3k* dE 

where S* = Si (k*) = and p = = 1.4500'lli with v* the neutron velocity at the 

critical energy for UCN production, E*, and c* the phonon group velocity at the critical 
energy. We have replaced JpS [q,uj = dcu where the subscript means the integration 
is carried out across the phonon peak with i3Si{k*) where 



Si{k*) = J S{k\uj)du (11) 

is the contribution to the structure factor from the one phonon scatteringi'El and (3 represents 
the increase in the integral due to the fact that the two curves (Q) and (|^) cross at an angle, 
so the path of integration is longer than in the calculation of (|llD, as was first pointed out 
by Pendleburylii. Numerical integration of the measured S [q,u! = ^g^)i across the phonon 
peak, (fig. 1) confirms this procedure. (Note that the factor of 2 in equation 9 of 0, which 
holds in the case of a linear dispersion relation should be replaced by /? = 1.45 for the real, 
non-linear dispersion relation) 
Thus (see equation |1^) 



Pf;civ = 9. 44 X 10-9^^^ UCN/ cm'/ sec (12) 
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-G: — Numerical calculaLiuiis uf the iiiulLiphuiiuii duwii scaLLeriiig cross secLiuii 



2 

Fig. |1] sliows S{q,uj = ^) as a function of q obtained from 0. We have extrapolated 
the data above 1.2 in fig. 0, taking into account the known value of 5* (g) for the 

one phonon and multiphonon contributionsiS. Measurement s0 show that S{ki,uj = ^) 
is essentially zero for ki >2 A"^. 

As an illustration we apply the above technique to various measured and calculated 
slow neutron spectra that are under discussion as locations for various superthermal UCN 
source projects. We will take kc= l-l x 10^^ A^^ corresponding to the critical energy for 
Beryllium. 
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1. North Carolina State proposed UCN source 



North Carolina State University researchers have proposed to locate a UCN source in 
the thermal column of the campus 1-MW PULSTAR reactor after removing the graphite. 
The expected average cold neutron flux in the UCN converter was estimated using MCNP 
simulations with a detailed model of the PULSTAR core and a conceptual model of the UCN 
source. The core's 4 control blades and 25 fuel assemblies, each consisting of a 5 x 5 array of 
fuel pins containing 4 % enriched U02, were faithfully represented. The model of the UCN 
source consisted of a cold neutron source surrounding a UCN converter about 10-cm deep 
in a tank of D20 in the thermal column void. The cold source moderator was 1-cm-thick 
cup-shaped solid methane (408 g, 22 K) contained in a 2-mm-thick-wall aluminum chamber 
(1164 g). The UCN converter was modeled by a disk of liquid ortho-deuterium (4-cm thick x 
14-cm diameter, 20 K) in a 5-mm-thick-wall aluminum container (484 g). Neutrons leaving 
a bare face of the reactor core were channeled into the D2O tank by a 44-cm diameter x 
75-cm long void in a beryllium assembly located between the reactor core and the reactor 
tank wall. The simulated spectrum is shown in fig. |I], along with the spectrum of some 

2 

other sources discussed below and S{q,uj = ^). 
For this spectrum 
/ $ (El) S{h,uj = ^)dki] = 5.22 X 10^ 
Substituting into equation ( 
(A^ = 2.18 X 10^2) 
we have 



P^p = N X 4nb'^ X 5.22 X 10^ X = 984 UCN/cm^/ sec 



1.3 
3 



For the one phonon contribution we have, taking ^ 

m 



from fig. ID and using equation 



ip 



9.44 X 10" 



dE„ 



9. 44 X 10"^ X 7.9 X 10^° = 745 UCN/cm^/ sec 



These figures will be reduced by a factor of 2 when the fiux (calculated here for a modest 
source volume as described above) is averaged over a more reasonable 20 liter volume. The 
effects of construction materials in a realistic design have not yet been taken into account. 
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The advantages of such a design were discussed in |T8[. The main point is that the UCN 
production process is independent of the direction of the incident neutrons so that use is 
made of the full solid angle of the flux. Installations at a beam position, as discussed in the 
next section, lose a solid angle factor of about 10^ [0 in comparison to installations inside 
a 4:71 flux, although such sources outperform any other type of UCN source that could be 
installed at the same beam position. 



2. Proposed UCN source at the Spallation Neutron Source (SNS) 

It has been proposed to place a superthermal source of UCN on a monochromatic 8.9A 
beam at a guide tube at the SNS.i It is interesting to ask what the loss of UCN production 
rate associated with the use of a monochromatic beam would be in the case of the simulated 
spectra available for the planned guide tubes. Spectra have been simulated for two types 
of guide: an 'ordinary' supermirror guide and a so-called 'ballistic' guide using tapered 
sections to focus and defocus the beam at the borders of a long free-flight region. This is an 
interesting case for the present discussion as the 'ballistic' guide has more flux at the critical 
8.9A wavelength but less flux at shorter wavelengths (A < 6A). The predicted fluxes for the 
two types of guide are shown in fig. |^, while fig. |I| shows a comparison of the SNS ballistic 
guide with the other sources considered in this paper. 

The result for the 'ordinary' guide is that 

.2 



SO that 



j $ (El) S{h,u= ^)dki = 5.38 X 10*^ 



Pmp = N X A-nh'^ X 5.38 x 10^ x = 1. 015 UCN/cm^/ sec 



For the one phonon production rate we find using 

Pip = 9. 44 X 10-^ 
For the 'ballistic' guide we have 



= 1.8 UCN/cm^/ sec 

E* 



J $ (^i) S(A;i, u = '^)dki = 5 X 10^ 



so that Pmp = 0.94 and the one phonon production rate is: 2. 36 UCN j err? j sec. 
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FIG. 2: Comparison of simulated spectra from the two types of guide proposed for the SNS 
fundamental physics position function of q, A ^ 

Thus we see that the losses associated with rejecting the higher energy neutrons in the 
beam by the monochromator are less significant than the losses (50%)^ connected with the 
monochromator and associated beam bender. 

3. HMI beam spectrum 

This spectrum was measured at the NL4 neutron guide at the Hahn Meitner Institut, 
Berlin by time of fiight. In this context we should remind the reader that the measurement of 
the spectrum from a neutron guide is quite complex as the spectra in general are very strong 
functions of the fiight direction and position in the beam cross section. For the measured 
spectrum (normalized to 1 at the peak of dip/dX) we find: 



j ^{E,)S{h,u = ^)dh 



.025 



yielding 



Pmp = 4. 72 X 10-9 



109 



1 .3 
1 .2 
1 .1 

1.0 H 

0.9 
0.8 
0.7 
0.6 H 
0.5 - 

0.4 - 

FIG. 3: 
0.3 - 

0.1 - 
0.0 - 



















— >K — HMI spectrum, normalised to max 
— • — UCN production rate, arbitrary units 
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Fig. ^ shows the spectrum and calculated UCN production rate as a function of the 
wavelength of the incident neutrons, A, similar to what would be observed in a time of flight 
measurement. 



II. DISCUSSION 

The results are summarized in the table below. The column labelled 'Maxwell" refers to 
a Maxwellian spectrum corrected for an ideal guide transmission r oc and cut off at 3.8A 
minimum wavelength. 

We see that the multiphonon contribution as expected is a rather strong function of the 
source spectrum varying from less than to slightly more than the one phonon production 
rate for the realistic spectra considered here. 

Table 1 - Predicted production rates 
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NC State* 


SNS ord* 


SNS ball* 


HMI a.u. 


Maxwell 


multi-ph 


490 


1.0 


.94 


4.7 


1.7 


single ph 


375 


1.8 


2.4 


5.5 


1.5 


mph/lph 


1.4 


.55 


.4 


.85 


1.13 



*UCN/cmVsec 



(The NC State figures have been reduced by the factor of 2 for a 20 liter source volume 
as discussed above.) 

Thus we see that the inclusion of the multiphonon production amounts to at most a 
little more than a factor of 2 increase in UCN production. In the case of a cold beam 
the multiphonon contribution is a small correction to the single phonon production so that 
the use of a monochromatic beam, which offers significant operating advantaged, would be 
accompanied by a minor loss in UCN production. 

In contrast, for sources where the Helium is exposed to the total thermal flux as in 
the North Carolina State proposal or at a dedicated spallation source, the multiphonon 
contribution can amount to slightly more than a factor of 2 increase in UCN production. 
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